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ABSTRACT 

We present a semi-automated method to search for strong galaxy-galaxy lenses in 
optical imaging surveys. Our search technique constrains the shape of strongly lensed 
galaxies (or arcs) in a multi-parameter space, which includes the third order (octopole) 
moments of objects. This method is applied to the Deep Lens Survey (DLS), a deep 
ground based weak lensing survey imaging to R ~ 26. The parameter space of arcs 
in the DLS is simulated using real galaxies extracted from deep HST fields in order 
to more accurately reproduce the properties of arcs. Arcs are detected in the DLS 
using a pixel thresholding method and candidate arcs are selected within this multi- 
parameter space. Examples of strong galaxy-galaxy lens candidates discovered in the 
DLS F2 field (4 square degrees) are presented. 
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1 INTRODUCTION 

Strong gravitational lensing is a powerful method with which 
to study the mass and mass distribution of foreground 
lens galaxies, independent of dynamical assumptions. The 
majority of co nfirmed galaxy sca le lenses are still galaxy- 
quasar lenses (|Falco et al.l 1 19991 ). though th e number of 
galax y-galaxy lenses is quickly catching up fi.e. lBolton et al] 
120071 ). Recent work in gravitational lens modeling suggests 
that galaxy-galaxy lens systems, where the source is ex- 
tended (a galaxy), p rovide the most accurate lens models 
ijWarren &: Dve|[2003l ). Galaxy-galaxy lens systems can also 
be used to study the lensed galaxies th emselves, providin g 
a window into higher redshift galaxies (|Allam et alj|2007l ). 
With futur e wide field surveys such as th e Dark Energy Sur - 
vey (DES) (|Annis et al.ll2005l ) or LSST l|Tvson et al.ll2002D . 
large samples of strong galaxy-galaxy lenses could also be 
used to put constrai nts on cosmolog y through th e statis- 
tics o f strong lensing (|Kochaneklll99rj ) (|Chaell2002l ) (|Linderl 
|2004) . 

Miralda-Escude & Lehar (1992) were the first to point 
out that a large number of Einstein rings (galaxy scale lens 
systems in which the source is a galaxy) should exist in the 
optical. They pointed out that the detection of rings depends 
on (1) angular resolution (2) the ability to identify rings 
among the large number of optical sources. They suggested 
looking for likely lenses (massive ellipticals) and subtracting 
out the light from the foreground lens galaxy, as a method of 
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detecting candidate ring systems. Many systems have since 
been detected in the optical, but mostly through visual selec- 
tion. Examples include the four systems in the HST Medium 

1 — 11 " 1 

deep survey (Ratnatunga et al. 1999), one system in the Ul- 
tra Deep Field (IBlakeslee et al.ll2004). and th e three systems 
in the AEGIS survey ilMoustakas et. aj|2007h . The GOODS 
survey (jFassnacht et all 20041 ) has discovered a number of 
candidate systems using a method similar to that suggested 
by Miralda-Escude & Lehar (1992). Here image subtraction 
in GOODS works well because of the HST resolution in all 
filters and relatively small survey area. Much progress has 
been made recently in finding strong galaxy-galaxy lense d 
systems in the SDSS redshift survey ( Bolton et al.ll2004l ). 
Here the emission line(s) of lensed source galaxies are de- 
tected in the spectrum of foreground lens galaxies. This 
search method is quite powerful but requires galaxy spectra, 
and therefore is not applicable to purely imaging surveys. 

Large area ground based optical surveys, especially 
weak lensing surveys, should be excellent places to search 
for candidate strong galaxy-galaxy lens systems. These sur- 
veys are wide area, deep, and because of the observational 
requirements to measure the weak lensing signal, have rel- 
atively good resolution (seeing). Good seeing is essential in 
separating the arc from the host lens galaxy, and the wide 
area increases the number of possible systems. Since weak 
lensing surveys are also focused on the d etection of galaxy 
clust ers, from maps of weak lens ing shear ([Kaiser fc Squires! 
Il995t ) or optical cluster finding l|Hansen et al.ll2005l )." struc- 
ture within the survey fields is well understood. Thus the 
environments of strong lens systems, which can potentially 
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affect the modeling of strong lenses dKeeton fc Zabludofj 
12004 ). can be taken into account. 

Semi-automated algorithms to search for strong galaxy- 
galaxy lenses are beginning to emerge, for instance 
in the Can ada-France Hawaii T elescope Legacy Survey 
(CFHTLS) (|Cabanac et al l 120071 1. Other arc finding al- 
gorithms exist but are currently limit ed to searching for 
cluster arcs at space based resolution (|Lenzen et al.ll2004l ) 
l|Seidel fc Bartelmannl 120071 ). The development of auto- 
mated algorithms to search for strong galaxy-galaxy lenses 
is of interest in order to create large uniform samples 
in future imaging surveys wh ere selection effects are un- 
derstood (|Cabanac et al.ll2007h . The current generation of 
weak lensing surveys such as the Deep Lens Survey (DLS) 
jWittman et alj|2002h and the CFHTLS can be used as test- 
ing grounds for the development of these search techniques. 

Here we present a new method to search for strong 
galaxy-galaxy lens systems in optical imaging surveys. Our 
technique focuses on using the shape of the strongly lensed 
background galaxy (the object which visually stands out in 
an image) to search for systems. This method has the ad- 
vantage of being independent of lens galaxy type, allowing 
us to recover systems that are due to massive ellipticals or 
potentially dark haloes. Systems with ellipticals as a lens 
galaxy are the most massive and will produce the highest 
separation images, making them easier to detect in ground 
based imaging. Systems produced by dark haloes should be 
detectable from the ground, if they do exist. Since the lens 
is dark in this case, there are no lens-arc separation issues 
to deal with. Excellent dark candidate systems would have 
a primary arc and secondary arc (with the same color as 
the primary) at a smaller separation from the lens centre. 
In general spiral systems will be harder to detect in ground 
based imaging since they are lower mass and therefore pro- 
duce images with lower lens-arc separation. We apply our 
technique to the Deep Lens Survey F2 field. 

The outline of the paper is as follows: In 321 we discuss 
the dataset used to apply our search technique. In Sj3] the 
properties of arcs are discussed, which are used to define an 
arc parameter space. Arcs are detected using a pixel thresh- 
olding method described in SJU In Sj5]the parameter space of 
arcs is generated using simulations of strong galaxy-galaxy 
lensing. Results of these simulations are discussed in 3H] In 
^7] we apply our search method to the DLS F2 field and 
present examples of detected candidates. 

Throughout the paper we assume a standard cosmology 
with Q m = 0.3 and Qa = 0.7. 



2 DATA 

Data us ed for this study are obtained from the Deep Lens 
Survey (jWittman et all 120021) a wide area (~ 20 degrees) 
multi-band (BVRz 1 ) weak lensing survey im aging to R ~ 
26. T he survey uses the MOSAIC imagers l|Muller et. all 
1998) on the Mayall 4-m at Kitt Peak and the Blanco 
4-m at Cerro Tololo, Chile. The survey is spilt into five 
separate 2° x 2° fields, with each field spilt into a 3 x 3 
grid of 40 arcmin x 40 arcmin subfields. The R band co- 
added images total 18, 000 seconds of exposure and the BVz' 
band total 12, 000 seconds. Details of the basic data reduc- 




Figure 1. The coordinate system in which the arcness of each 
object is measured. The image coordinates are given by the (x,y) 
axes, and (a,b) define the object major and minor axis. The po- 
sition angle (8) of the object is measured relative to the image x 
axis. The arcness of the object is chosen to be measured relative 
to the major axis of the object. 



tion, calibration, and co-addition pipelines are described in 
IWittman et~al1 (|2006l ). 

In this study we restrict our search to the DLS F2 field 
(centred on a = 09 h 18 m 00 s 5 = 30°00'00" J2000) which has 
complete imaging. For our strong lens search we use only 
the R band since this is the deepest filter with the best see- 
ing (< 0.9 arcsec). The four color filters in the DLS will be 
used in future weak lensing studies to derive photometric 
redshifts for source galaxies. At the time of our study pho- 
tometric redshifts were not available for the DLS F2 field, 
and therefore are not used in our analysis. 



3 SEARCH TECHNIQUE 

The primary motivation for our search technique was to 
measure what visually stands out to the eye when look- 
ing at a strong lensed galaxy, that is the arc shape. It was 
first pointed out by Goldberg & Natarajan (2002) that the 
shape of an arc is stored in one component of the octopole 
(or third order) moment. Th eir study and subsequent stud- 
ies l|Goldberg fc Baconll2005l ) focused on examining the oc- 
topole moments in the weak lensing limit. Here we use this 
component (along with the other third order moments) to 
search for arcs in the strong lensing regime. 



3.1 Octopole moments 

Since an arc can have a random orientation in an image, 
the octopole moments are measured relative to each object 
instead of the image axis. For each object we measure the 
octopole moments relative to the object major axis (Fig.[l|. 
By definition the position angle (ff) of the object is the angle 
between the x image axis and the major axis of an object, 
where 




Figure 2. Examples of simulated strong galaxy-galaxy lens sys- 
tems with different measured levels of arcness. In each, the data 
image is mapped to the segmentation image at the threshold in 
which the arc is detected. The foreground lens galaxy is centred 
on each image, and the arcness is measured on the primary arc 
in each system. From left to right the arcness of the primary arc 
in each system is : 0.20, 0.15, 0.10, 0.05. 
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3.2 Ellipticity 

In addition to the octopole moments, the ellipticity of a 
strongly lensed galaxy is also a useful arc property. The el- 
lipticity is measured in terms of its ei and ei components, 
where ei = * x ~ vv e2 — 21 *y — anc i the total elliptic- 

1 xx\ !yy *xx\lyy 

ity (e) is given by e = -y/ef + ef . The quadrapole moments 
are measured here in the most basic sense, that is I xx — 

^Zih^i-*) 2 r _ T.ih(Vi-v) 2 _„ j r _ S, h(x t --)(y,-y) 
J*. I f ' W ~ Y, • !% ' X V — Y !i 

No weights are used in calculating the second order moments 
in order to avoid making assumptions about the light profile 
along the arc. 
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The octopole moments are measured in the frame of each 
object using the rotation matrix 
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where x and y define the centroid of the object, and 9 
is the position angle. The four components of the oc- 
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topole moments are then, l a ab = y j - — ' = 
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where a and b are the semi-major and minor axis. The 
octopole moments are further made unitless and size in- 
dependent by normalizing by the cube of the object size 
(size= \JI XX + lyy)- The quantity I xx + I yy is rotationally 
invariant but we have measured this in image coordinates. 

Relative to the major axis, I aa b is the component of the 
octopole moment that measures the arc shape of the strongly 
lensed galaxy. To provide a concavity independent measure 
of the arc shape we take its absolute value |7 a at|- This nor- 
malized component is referred to throughout the rest of the 
paper as the galaxy "arcness" (given by the symbol : a), 
where 

\Iaab\ 
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To better illustrate the arcness of strongly lensed galaxies, 
examples of simulations with different values of measured 
arcness are shown in Fig. [2] Arcs with a high level of arcness 
correspond to systems in which the system is significantly 
tangentially distorted. 

Strongly lensed arcs measured relative to the major axis 
should also have a small I a t,t value, making this component 
of the octopole moment useful as a control. This component 
is referred to as the "anti-arcness" (given by the symbol : 
o), where 
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The remaining components of the octopole moments 
Iaaa and Ibbb, measure the skewness along the major (s a ) 
and minor (st) axes respectively 
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3.3 Arc S/N 

Another important property of arcs to consider is their range 
of expected S/N. The object S/N is computed from the stan- 
dard S/N equation using the object flux, the noise per pixel, 
and the gain (~ 54e~pixel _1 , for the co-added DLS R band 
image) at each detection threshold. At some minimum S/N 
it is expected that the measured arcness becomes unreliable. 
The minimum S/N for the DLS is discussed in £16.31 



3.4 Additional parameters 

In addition to the above parameters we also tested others 
that were not particularly useful. These included the arc 
length and arc width, which are typically used to study the 
shapes of cluster arcs jMiralda-Es cudej] 1 1993f) . At galaxy 
scales, arcs are not as drastically elongated as arcs produced 
by clusters, and most of the information contained in these 
parameters is stored in the ellipticity parameter. 



4 ARC DETECTION 
4.1 Pixel thresholding 

Detecting strongly lensed galaxies (arcs) is challenging since 
arcs typically occur near a foreground lens galaxy and also 
can have a complicated surface brightness profile. Their de- 
tection therefore depends on the ability of a detection algo- 
rithm to separate the arc from the lens, and also to avoid 
splitting the arc into multiple objects (and therefore lose 
its arc shap e). We chose to detect arcs in the DLS using 
SExtractor l)Bertin fc Arnou"tslll996t ). a object detection and 
shape measurement program. SExtractor typically detects 
objects by splitting a group of connected pixels using a 'de- 
blending' algorithm which assigns pixels to each object by 
computing a probability of belonging to each object. In the 
case of separating a lens galaxy and an arc, this step can 
possibly affect the shape of the resulting arc. Since there is 
no way for the user to control this step, we instead chose 
to use a pixel thresholding approach. Here the deblending 
algorithm in SExtractor is turned off and the program is run 
at a series of detection thresholds. This is similar to how the 
SExtractor deblending works, except here the shape infor- 
mation is preserved at each threshold, and only thresholding 
determines which object a given pixel is assigned. 
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4.2 Implementing pixel thresholding 

For a given detection threshold, SExtractor is used to cre- 
ate a segmentation image where the detection threshold is 
measured relative to the image RMS value. The segmenta- 
tion image is then used as input to our custom cataloging 
software which measures the octopole moments. Using this 
global per pixel value of the noise is better than using the 
SExtractor background map, which is a smoothed map of 
the image noise. We additionally made minor modifications 
to version 2.3.2 of the SExtractor code in order to accept a 
large number of objects in the segmentation image, which 
primarily occurs in the DLS fields at low detection thresh- 
olds. 



4.3 Thresholds 

The minimum object area is set sufficiently small (below the 
PSF) in order to avoid rejecting real arcs in the DLS. The 
minimum detection threshold to use was determined from 
fOO realizations of the simulation noise, where the threshold 
was increased until there was a zero probability of detecting 
noise as a real object. This resulted in a minimum detection 
threshold of 2.5. The maximum detection threshold is not 
as well defined as the minimum, but we set this detection 
threshold relative to the maximum S/N arc produced by 
our simulations, which is separated from the lens galaxy at 
a detection threshold of 20. 

Logarithmically spaced thresholds are used from the 
minimum detection threshold up to the maximum thresh- 
old. As the number of thresholds is increased we expect to 
reach a value where there is no gain in using more thresh- 
olds. For a large set of simulations we show in Fig. [3] the 
fraction of detected simulated arcs with arcness a > O.fO as 
a function of the number of detection thresholds. The frac- 
tion of detections here flattens for ~ f 5 detection thresholds. 
This is the number of thresholds at which the simulations 
and data are analyzed. The overall fraction here is low since 
our simulations produce systems with all values of arcness, 
which are dominated by systems with a low arcness level. 



5 STRONG LENS SIMULATIONS 

To estimate the parameter space of arcs in the DLS we use 
a large number of simulations (90, 000) of strong galaxy- 
galaxy lensing. This is a sufficiently large number of sim- 
ulations to estimate the arc parameter space in the DLS 
while still being computationally feasible. To create realistic 
properties for the arcs in each simulation we u se real galax- 
ies fr om deep HST fields for our sources (as in lHoresh et al.l 
120051 ). instead of artifical sources. 

5.1 Source galaxies 

5.1.1 Source image database 

To create realistic shapes for strongly lensed galaxies (or 
arcs) real galaxies from deep HST exposures are used for 
our sources. The deepest HST exposures from which to 
select sou rce galaxies were the Hubble Deep Field North 
(HDFN) (|Williams et all If996l ). the Hubble Deep Field 
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Figure 3. The fraction of simulated arcs detected above an ar- 
cness level a > 0.10, as a function of the number of detection 
thresholds. As the number of detection thresholds increases, the 
fraction of recovered simulated arcs flattens at ~ 15 thresholds. 
This is the number of detection thresholds used to sample the 
shape of arcs in the DLS. 

South (HDFS) l|Casertano et alJl200Ch . t he Great Obser- 
vatories Origins Deep Survey (GOODS) llGiavalisco et al.l 
| 2004) . and the Ultra Deep Field (UDF) jBeckwith et al. 
2006). From each field, galaxies which have the highest S/N 
were picked from the F606W filter (which matches reason- 
ably well with the DLS R band). Each galaxy was thresh- 
olded at 4<r (per pixel) and cleaned to create a "noise free" 
image. This resulted in a database of 90 unique galaxies, 52 
from the UDF, 19 from GOODS, 9 from the HDFN, and 
10 from the HDFS. Examples of these galaxies are shown in 
Fig. [4] At the time we began our study morphology was not 
available for galaxies from the UDF, so galaxies were visu- 
ally classified into three categories: spirals, ellipticals, and 
irregulars. The HDF's did have morphology available, but 
to be consistent with the UDF sample, Sbc's and Scd's were 
grouped together as "spirals". For each run of the simula- 
tion, a source galaxy in the database is selected and rotated 
to a position angle (both at random), to reflect the random 
orientation of source galaxies. The source galaxies in the 
database were also photometrically scaled to the simulated 
lens galaxy zeropoint. 

5.1.2 Source image reseating 

In using real galaxies for sources we are limited by the 
total number of galaxies with sufficiently high S/N avail- 
able from the HST fields. The effective number of available 
source galaxies is increased by rescaling each initial source 
galaxy image to reflect the size, magnitude, and redshift dis- 
tribution as determined from HST surveys. To rescale the 
images we initially select a new magnitude for the source 
galaxy at random from a power law distribution in the range 
22 < m sourco < 25. The flux in the initial database source 
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Figure 4. Examples of source galaxies in the database used in 
our simulations of strong galaxy-galaxy lensing (Left: Spiral, Mid- 
dle: Irregular, Right: Elliptical). Galaxies were selected from the 
F606W band of the Hubble Deep Fields (North and South), the 
Ultra Deep Field, and the Great Observatories Origins Deep Sur- 
vey. Each galaxy was thresholded at 4ct and cleaned to create a 
"noise free" version of the galaxy. 



galaxy image is then scaled to the new source galaxy mag- 
nitude. A new size is assigned to the source galaxy based on 
the magnitude-size distribution of galaxies from the UDF. 
From the new magnitude and the associated magnitude bin 
in the UDF catalog, a new size is randomly assigned to 
the galaxy and the image is scaled to the new size using 
a linear interpolation routine. We filtered the initial UDF 
catalog to include only galaxies where the size is well mea- 
sured (R 1/2 > 1.5 x PSF, where the PSF ~ 0.06 arcsec 
and R1/2 is the half-light radius). A redshift is assigned 
the source galaxy using the magnitude-redshift distribution 
of galaxies from t he combined photometric re dshift cata- 
logs of the HDFN {Fernandez-Soto et al.|[l999l) and HDFS 
ijYahata et alJ [2000). The source galaxy is randomly as- 
signed a new redshift based on its initial morphological type 
and the associated magnitude bin in the redshift catalog. 
The initial redshift catalog is separated by galaxy morpho- 
logical type, and filtered to include galaxies within the range 
0.2 < z < 3.0. We chose a maximum source redshift cutoff of 
z = 3.0 since few galaxies are expected to be detected above 



this redshift in the DLS R band at this magnitude and size 
limit. 



5.2 Lens galaxies 

Since the lensi ng cross section is expecte d to be dominated 
by ellipticals l|Fukugita fc Turnen Il99 ll ) , we use elliptical 
galaxies for our lens galaxies. For each simulation a lens 
galaxy magnitude is randomly chosen from a power law dis- 
tribution in the magnitude range 21 < mi cns < 23. The 
lens galaxy is assigned a redshift based on the magnitude- 
redshift relation of elliptical galaxies determined from HST 
fields (J. A. Tyson 2007, private communication). Photomet- 
ric redshifts of lenses in the DLS are not used since these 
were not available during the course of this study. Knowl- 
edge of the exact lens redshift distribution in the DLS is 
not crucial since this has only a small effect on the angular 
diameter distance ratio in equation J5). 

For each lens galaxy a velocity d ispersion is esti- 
mate d using the Faber- Jackson relation (|Faber fc Jackson] 
1 19761 ). We consider only galaxies with a velocity disper- 
sion 250kms~ 1 < a v < 400kms _1 , since this is the range 
which will produce systems with measurable image separa- 
tions (> 1.0 arcsec) in the DLS. To simulate the effects of 
separating the arc and lens galaxy in real DLS data, an el- 
liptical lens galaxy is included in the final simulated image. 



5.3 Gravitational lens mapping 

Each source galaxy image is gravitationally lensed using ray 
tracing via the lens equation 6 = /3 + a, where is the 
image position, f3 is the source position, and a is the deflec- 
tion given for the projected scalar potential. The deflection 
angle a. is related to the scalar potential via a — Vi/>. To 
increase the resolution of the lens mapping, each pixel in 
the lens plane is subdivided into a 4 x 4 grid of sub-pixels. 
For the lens pote ntial we chose to use a project ed scalar el- 
liptical potential IjBlandford fc Kochaneklll987l ) because of 
its analytic simplicity and since it generates realistic lensed 
images. The deflection for each sub-pixel is calculated from 
the x and y components of the potential given by 



ip(x,y) = i>o\/l- 
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Here Di 3 is the angular diameter distance between the lens 
and source, D s is the angular diameter distance of the 
source, and a v is the velocity dispersion of the lens. The 
position of the centre of the lens is given by a;i ens and yiena, 
and r c is the core radius. Since we are only interested in the 
primary arc in each system (which is sensitive to the total 
enclosed mass) we have chosen a small core radius for each 
system. The lens galaxy ellipticity (eiens) is coupled to the 
potential ellipticity (e p ) via the function 



0.2e,, 



0.6 + eie 



(11) 



which constrains the potential ellipticity to the range < 
e p < 0.2, where the surface density remains physical. Each 
system is ensured to be strongly lensed by placing the source 
within the Einstein angle for the system. An example of 
lensed source galaxy is shown in Fig. [5] (Right). 



5.4 Scaling to the DLS 

Each strongly lensed image (0.03 arcsec) is scaled to the DLS 
pixel scale (0.257 arcsec) using a linear interpolation routine, 
where the flux is conserved in the scaling. The lensed source 
galaxy image is added to the lens galaxy image to create the 
strongly lensed system. An example of a lensed system scaled 
to the DLS pixel scale is shown in Fig. [6] (Left). After scaling 
to the DLS pixel scale, the image is set to the DLS seeing 
(0.9 arcsec) by convolving the image with a Gaussian kernel, 
where the sigma of the Gaussian is related to the FWHM 



by a = 



An example of a lens system smoothed to 



2^/2111(2) ' 

the DLS seeing is shown in Fig. [6] (Right). 

To set the S/N of the simulations to the DLS S/N, an 
object in the simulation (after convolution) was matched to 
an object in the DLS with the same magnitude and size. 
For the matching object in the DLS, the noise (odls) was 
measured in the region surrounding the object. The noise 
to set in the simulation {a S vm) is determined by <j S i m — 



6 Kubo & Dell 'Antonio 




Figure 5. {Left): A galaxy selected from the simulation database is scaled to a new size, magnitude, redshift, and shifted to an impact 
parameter. (Right): The simulated galaxy on the left is gravitationally lensed by an elliptical potential. The structure within the source 
galaxy is evident in both images each of which have a pixel scale of 0.03 arcsec pixel -1 . 



Figure 6. (Left): The image of the gravitationally lensed source galaxy and foreground lens galaxy is scaled to the DLS pixel scale 0.257 
arcsec pixel - and a foreground lens galaxy is added. (Right): The system on the left is convolved with the DLS PSF using a Gaussian 
with FWHM = 0.9 arcsec. 



tdls f^ls ' wnere -fsim and Fdls are the flux of the object 
in the simulation and data respectively. 

For each simulation Poisson noise is generated on a 
background (B) which is set by 

p _ (cr s im x gain) 2 



where a new random seed is used for each simulation. For 
the gain we use the same value as in £13.31 The background 
value is subtracted from the noise image to create the final 
noise image. A simulated strong galaxy-galaxy lens system 
with noise is shown in Fig. [7] 



6 SIMULATION RESULTS 
6.1 Detecting simulated arcs 

We use an automated method of extracting arcs from the 
simulations since a large number of simulations are gener- 
ated (~ 90000). Arcs are detected in the simulations at ex- 
actly the same thresholds that are used for detection on the 
DLS data. Since the thresholded catalogs for each simula- 
tion contain a galaxy lens and at least one arc, a minimum 
of two objects is required to be detected in each catalog. 
Objects are searched for which have a minimum separation 
of 1.0 arcsec from the centre of the lens and are resolved 
(V2xsize > 1.00 arcsec). Detected objects meeting these cri- 
teria in each thresholded simulation catalog are considered 



Figure 7. A completed simulation of a strong galaxy-galaxy lens 
system. In the image noise has been added to the combined image 
of the foreground lens galaxy and gravitationally lensed source 
galaxy. The noise scales the system to the correct DLS S/N. Most 
of the structure within the DLS scaled arc is gone, except for the 
central region of the source galaxy. 

arcs. For a given simulation, if an arc is detected at multiple 
thresholds one of the thresholds is randomly picked as the 
shape of the arc and used for the result of the simulation. 

6.2 Arc parameter distributions 

The detection method described in £|6 . 1 1 was used to de- 
termine the properties of our simulated arcs. In Fig. [8] 
the normalized distribution of simulated arcness is shown 
(solid line), where the peak in the distribution occurs for 
a small value of arcness. The low end of arcness represents 
source galaxies that have been strongly lensed, but where 
the galaxy shape is not particularly curved. Data for a sin- 
gle subfield is also shown (dot-dash line) where the peak in 
arcness for data occurs at a lower value of arcness. The off- 
set between the simulations and data indicates that there is 
discrimination power in the arcness parameter, between arcs 
and other objects in the subfield. The dashed line indicated 
at arcness value of 0.10 represents a possible cut in arcness 
above which to select arcs which are significantly distorted. 

In Fig. [5] the distribution of anti-arcness is shown for 
arcs with an a > 0.10. These highly distorted arcs are dis- 
tributed tightly around zero in this parameter space. This is 
because arcs which are significantly distorted relative to the 
major axis are not significantly distorted relative to their 
minor axis. The resulting distributions for the remaining 
octopole moments, the skewness along the major and mi- 
nor axes, are similarly distributed about a zero value. Both 
distributions are Gaussian with a FWHM~ 0.8 for s a and 
FWHM~ 0.03 for Sb, and therefore the skewness along the 
major axis varies significantly more than the skewness along 
the minor axis. 

In Fig. \W\ the distribution of ellipticity is shown for 
simulated arcs with a > 0.10. The majority of significantly 
distorted arcs occur at high ellipticity e > 0.6, where a line 
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Figure 8. The normalized distribution of arcness for simulations 
(solid) and for data in a DLS subfield (dot-dash). Both histograms 
show objects which are extended (e > 0.6) and are above the 
minimum S/N level (S/N > 40). Most extended objects in a DLS 
subfield have a small arcness with the highest fraction of arcs at 
a ~ 0.02. The most probable arcness for simulated arcs occurs 
at a higher arcness level a ~ 0.06. The dashed line at a = 0.1 
represents a possible cut to make in the simulations above which 
to select arcs which are significantly distorted. 

0.4 1 1 1 1 ' ' 1 1 1 1 ' 1 1 1— 1 
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Figure 9. This normalized distribution of anti-arcness (a) for 
simulations which have been significantly distorted (a > 0.10). 
The anti-arcness values are peaked around zero, which indicates 
that arcs which are significantly distorted relative to their major 
axes have a small distortion relative to their minor axis. 
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Figure 10. The distribution of ellipticity for arcs which are sig- 
nificantly distorted (a > 0.10) and which meet the minimum S/N 
(S/N > 40). The most probable ellipticity occurs for arcs with 
(e ~ 0.75). The dashed line at an ellipticity of e = 0.6 indicates 
a cut to make in ellipticity, above which to select arcs that are 
extended. 



at this ellipticity is shown to represent a cut above which to 
select arcs that are extended. 



6.3 S/N 

In Fig. [TTJ (Left) the distribution in S/N for arcs with 
an a > 0.10 is shown. The most probable S/N occurs 
at S/N ~ 70 with virtually all simulated arcs having a 
S/N < 1000. The minimum simulated arc S/N extends to 
very small values of S/N, however we expect that at some 
S/N level the measured arcness becomes unreliable. To de- 
termine this minimum S/N, a set of simulated arcs were 
smoothed to seeing of 0.9 arcsec, and the arcness for each 
arc was measured at a series of different background levels. 
For this bootstrap method, 100 realizations of the noise were 
used at each background level. The results of this test are 
shown in Fig. [TT] (Right). To compare arcs from different 
simulations, the arcness at all S/N levels are scaled relative 
to the highest S/N level (which has been scaled to an arcness 
of 0.10). Measurements at different noise realizations for the 
arc are then grouped into different bins of S/N. Each point 
in the plot represents the median value for all simulations 
within that S/N bin. The error bars are relative to the arc- 
ness values at 20% (lower error bar) and 80% (upper error 
bar) within each bin. At high S/N levels the arcness for the 
arcs remains relatively stable, but the uncertainty in arcness 
becomes high for simulations with S/N 40. Therefore we 
set our minimum S/N to reliably measure the arcness at 
S/N = 40, indicated by a line in Fig. [TTJ (Left). We could 
have been more restrictive and raised the minimum S/N 
here, but did not want to reject a large number of potential 
strong lens systems with this cut. 



Table 1. Fraction of Simulations With Different Arcness Levels 



arcness range 


f (%) 


0.05 - 0.10 


15.4 


0.10 - 0.15 


10.9 


> 0.15 


6.7 



6.4 Lensing fraction 

Our simulations were generated by placing the source galaxy 
at random positions within the Einstein radius of the sys- 
tem. Probing the entire space of possible arc systems allows 
us to determine the fraction of systems that will produce a 
certain arc configuration. This fraction is shown in Table [TJ 
for the different ranges of arcness used in SjZl The fraction 
represents the number of simulations which are above the 
arcness level, extended (e > 0.6), and meet the minimum 
S/N level (S/N > 40). It is normalized by the total num- 
ber of simulations which are considered strongly lensed and 
meet our minimum lens-arc separation (> 1.00 arcsec). 

Arcs which are the most significantly distorted corre- 
spond to only a small fraction ~ 7% of the possible strongly 
lensed systems in the DLS; however as the arcness level is re- 
duced more of the parameter space of arcs is recovered. For 
example, including all arcs with a > 0.10 would correspond 
to ~ 18% of the possible parameter space, or including all 
arcs with a > 0.05 would correspond to 33% of the param- 
eter space. 



7 ARC SEARCH 

7.1 Appication to DLS data 

We searched for arcs in the DLS F2 field using the arc pa- 
rameters outlined in Sj3] independent of the position of any 
foreground lens galaxy. Each DLS F2 subfield is cataloged 
using the pixel thresholding method described in [J4] In each 
thresholded catalog definite spurious objects occur near sub- 
field edges and bright stellar haloes. These regions are fil- 
tered by creating SAOImage DS9 region files by hand. Using 
the probability selection described in H7.2I we performed a 
blind search for systems in the DLS F2 field ^7.31 Searching 
the entire parameter space of arcs produced a large number 
of spurious detections ~ 40, 000 per square degree. This is 
because all objects which are strongly lensed are searched 
for including those with a small shape distortion. To re- 
duce the number of spurious detections we restricted the 
type of arc to search for, searching only for galaxies which 
are significantly distorted and extended. This corresponds 
to searching for objects which have a high arcness and high 
ellipticity. The space of highly distorted and extended arcs is 
generated by selecting only these systems from simulations. 

Our blind search for high arcness systems produced a 
list of candidates systems in the DLS F2 which were visually 
inspected. Here objects which were obvious superpositions 
of two objects were rejected. We also rejected objects which 
were clearly spiral arms of a galaxy or were clearly associ- 
ated with another foreground object. Examples of typical 
spurious detections are shown in Fig. 1121 
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Figure 11. (Left): The S/N range of simulated galaxy-scale arcs in the DLS. The dotted line at a S/N of 40 indicates a minimum value 
of S/N that an arc should have in order for the arcness to be considered reliable. The most probable S/N occurs at S/N ~ 70 with 
virtually all simulated arcs having a S/N < 1000. (Right): A set of simulated arcs is set to different background levels and scaled to an 
arcness of 0.10. Each point is the median value within each S/N bin and error bars are set to the 20% and 80% values within each bin. 
The uncertainty in arcness becomes large at a S/N ~ 40 shown as a dotted line. 






Figure 12. Examples of typical spurious detections rejected during visual inspection. Centred on each image are objects which lie in the 
parameter space with a > 0.15 and e > 0.6. These types of objects are rejected upon visual inspection since they are clear superpositions 
of two galaxies. The dimensions of the above images are 26 arcsec X 26 arcsec. 



7.2 Probability selection 

Using the simulation results for each arc parameter in £|6 . 2 1 
and i]6.3l we calculated the probability that an object in the 
DLS is consistent with the properties of a strongly lensed 
arc. We desired a method which was general enough so that 
any new parameters could be added relatively easily. This 



ruled out a Bayesian approach because many of our arc pa- 
rameters are dependent on each other, and adding a new 
parameter would become cumbersome. The method chosen 
to calculate the probability was based on the occupation 
probability in the multi-parameter space defined by our sim- 
ulations. 

The number of arc parameters and the number of bins 
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Figure 13. Examples of arc candidates detected in the DLS F2 field. The scale of each image is 26 arcsec X 26 arcsec and the orientation 
is East is down, North is to the Left. The candidate on the (Left) DLS01 was detected in the parameter space a > 0.15 and e > 0.6. The 
candidate on the (Right) DLS02 was discovered in the parameter space 0.10 < a < 0.15 and e > 0.6. Both arc candidates appear centred 
towards a corresponding elliptical lens galaxy. 



for each parameter are used to define the multi-parameter 
space. The catalog of simulated arc shapes is used to ini- 
tially populate the space and each bin is normalized in or- 
der to assign a bin probability. The simulation probability 
bins are sorted in descending order while keeping track of 
the location of each bin. We additionally define a "proba- 
bility contour" for the space by summing from the maxi- 
mum probability until a cut in total probability is met. This 
probability contour can be used to reject large portions of 
the parameter space of arcs, however in general we keep the 
probability high in order to detect a large number of arcs. 
The bin sizes are made small in order to maximize the re- 
jection of spurious objects. To ensure that the bin sizes are 
not too small (and thus have a fragmented space) the grid 
in the multi-parameter space was required to be connected. 
A given parameter bin must be connected to the other sim- 
ulation bins in at least one direction along any parameter. 

7.3 Blind search 

In our blind search, arcs with high ellipticity were searched 
for at two different levels of arcness. The first level was for 
arcs with a > 0.15 and e > 0.6. Restricting the parame- 
ter space here reduced considerably the number of spurious 
detections, to an average of ~ 410 per square degree elim- 
inating 99.9% of all objects in a DLS subfield. One inter- 
esting candidate (DLS01) is discovered at this arcness level 
shown on the Left in Fig. [13] (North is to the left, East is 
down). The candidate arc has a measured arcness a = 0.18, 
ellipticity e = 0.73, and magnitude R = 23.5. The arc also 
appears centred relative to an elliptical lens galaxy located 
southwest of the arc with magnitude R = 20.37. 

We next searched for arcs in the parameter space with 
0.10 < a < 0.15 and e > 0.6. Within this range of arcness 



the number of spurious detections increases to ~ 1100 per 
square degree. Another interesting arc candidate (DLS02) 
is detected in this lower arcness space shown on the Right 
in Fig. [T3] (North is to the left, East is down). This arc 
candidate has a measured arcness a = 0.13, ellipticity e = 
0.80, and magnitude R — 23.9. The arc also appears centred 
on an elliptical lens galaxy located directly north of the arc 
with magnitude R = 21.96. 

Using our previously discussed criteria for an excellent 
candidate dark lens in SJT] we find no candidate dark lens 
systems in the DLS F2 for either parameter space. 

7.4 Visual search 

In addition to our blind search we also performed a visual 
search in the DLS F2. Here we viewed each subfield image 
at high contrast in order to look for arcs within bright lens 
galaxy haloes that our previous search could have missed. 
We recovered one arc candidate (DLS03) in F2 shown in Fig. 
1141 The system was missed in our blind search because pixel 
thresholding was not able to successfully separate the arc 
from the foreground lens galaxy at any threshold. This was 
due to the brightness of the foreground lens galaxy (R = 
18.19) and the extent of its halo in the DLS imaging. 

7.5 Candidate discussion 

Properties for each candidate strong lens system are given in 
Tabled The candidate system DLS01 has a large estimated 
Einstein radius (4.9 arcsec) which is likely due to its proxim- 
ity to nearby galaxy clusters. The preliminary weak lensing 
convergence map for the DLS F2 field (Wittman et. al 2006) 
reveals a large structure with three overlapping peaks. A red- 
shift survey of this field has identified these peaks to be Abell 
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Figure 14. The arc candidate DLS03 which is visually detected 
in the DLS F2 field. The scale of each image is 26 arcsec X 26 
arcsec and the orientation is East is down, North is to the Left. 
This system is not detected using our search technique because 
of the bright halo of the foreground lens galaxy (R = 18.19) in 
the DLS imaging. 



781, CXOU J09205 3+3029880, and CXOU J092110+302751 
ijGeller et. alll2005T ). Relative to the centre of the peak in the 
convergence map DLS01 lies 2.9 arcmin away. The visually 
detected candidate DLS03 also lies close to these clusters 
(7.4 arcmin from the centre) and has the largest estimated 
Einstein radius (6.7 arcsec) in our sample. The remaining 
candidate DLS02 has the smallest Einstein radius (3.4 arc- 
sec) and is separated by at least 20 arcmin from any peak 
in the convergence map. 

We have started follow-up of the visually detected sys- 
tem but leave a detailed discussion and lensing analysis for 
a future paper. Redshifts of the lens galaxies in the other 
systems are unknown, though future photometric redshifts 
in the DLS F2 will provide estimates. Spectroscopic follow- 
up will allow for definitive confirmation of these systems, 
however since the lens and arc in both of these systems are 
faint this will require follow-up on 8-m class telescopes. 

It is possible that the candidate arcs in DLS01 and 
DLS02 could be galaxies which are not strongly lensed but 
instead have a high intrinsic arcness, for instance due to a 
tidal interaction. Follow-up spectroscopy of our candidates 
would be needed to further examine this possibility. The fre- 
quency of objects with high intrinsic arcness in the DLS or 
any optical imaging survey is not understood, and warrants 
future study. 

Additional candidate systems with smaller lens-arc sep- 
aration could potentially be recovered by restricting our 
search to regions near foreground lens galaxies and subtract- 
ing out the lens galaxy light profile. For computational rea- 
sons we did not pursue this in our current study but this 
will be explored in future work. 



8 CONCLUSION 

A new method of searching for strong galaxy-galaxy lens 
systems in optical imaging surveys is presented. Our search 
method uses the third order moments of galaxies (along 
with the other arc parameters we have outlined) in order 
to identify candidate strong lens systems. Since detection 
depends only on the properties of the lensed source galaxy, 
this method can potentially detect systems independent of 
lens galaxy type, including systems produced by dark lenses. 
This method was applied to the Deep Lens Survey, a current 
wide field ground based optical survey. Within this dataset 
our method is sensitive to more massive systems where the 
lens-arc separation is > 1.0 arcsec. This method is partic- 
ularly successful on systems in which the foreground lens 
galaxy is faint, which potentially lie at higher redshift. Sys- 
tems near bright foreground lens galaxies (which likely lie at 
lower redshift) are harder to detect with this method, but 
are recoverable within a survey of this size by viewing the 
co-added images at high contrast. In the future other arc de- 
tection methods, instead of the pixel thresholding method 
used here, could help to improve the automated detection 
of arcs near bright lens galaxies. 

Our method should be useful for upcoming wide area 
ground based surveys such as the DES or future space 
based imagin g surveys such as SNAP (Supernova Acceler- 
ation Probe) (|Alderindl2005l) . Imaging in the DES will be 
shallower than the DLS, decreasing the number of spurious 
detection per square degree. Here the combination of wide 
area and good seeing should provide a dataset where a large 
number of strong galaxy-galaxy lens systems system could 
be detected. A significant number of systems could poten- 
tially be detected with SNAP, in particular small separation 
systems (< 1"), given the small space based point spread 
function. Current HST archival data would provide a good 
testing ground for our technique at spaced based resolution 
and should be explored. 
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Table 2. Candidate Strong Lens Galaxies in the DLS F2 Field. The Einstein radius is estimated 
from the R band data assuming the arc position traces the critical curve. Arcness and ellipticity 
parameters are not reported for the candidate DLS03 since this candidate was visually detected. 



Name 


R.A.(J2000) 


DEC.(J2000) 


■^4 ens 


R-arc 


6>E cst (arcsec) 


arcness 


ellipticity 


DLS01 


09:20:19.78 


+30:28:31.81 


20.37 


23.5 


4.9 


0.18 


0.73 


DLS02 


09:16:29.71 


+29:53:40.46 


21.96 


23.9 


3.4 


0.13 


0.80 


DLS03 


09:19:35.06 


+30:31:56.63 


18.19 


21.1 


6.7 


N/A 


N/A 
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